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Partial-wave analyses (PWAs) have a long history in the fields of vrA^ and AW scatter- 
ing. Due to the poor quality of low-energy antiproton beams and the resulting absence 
of accurate experimental data, analogous model-independent studies of the much more 
complex pp system have in the past always been impossible. In recent years, however, 
experimental progress has been very significant, in particular due to the coming in 1983 of 
the Low-Energy Antiproton Ring (LEAR) facility at CERN. While in the pre-LEAR era 
^ I spin-dependent observables and charge-exchange {pp —>■ nn) data were almost nonexistent, 
^ I the situation between 400 and 925 MeV/c is now quite good: the LEAR collaborations 
PS172, PS173, PS198, and PS199 have measured a variety of observables with impressive 
^ . accuracy. High-quality analyzing-power data have been obtained for the elastic [jl| and 
^ ! charge-exchange reactions. Very recently, even charge-exchange depolarization data 
have become available 0. Unfortunately, the practical difficulties involved in construct- 
ing a high-quality "cooled" antiproton beam of lower momentum are large. Consequently, 
the pp database below about 400 MeV/c is still by far not as good as one would like, in 
striking contrast to the pp case where very accurate data exist as low as = 0.35 MeV 
{Pl = 25 MeV/c). It also remains an outstanding experimental challenge to construct a 
polarized antiproton beam to further probe the spin structure of the interaction. 

During the last 10 years a new method has been developed by the Nijmegen group 
to perform PWAs of the abundant and accurate NN {pp and np) scattering data below 
Tl = 350 MeV g |]. With the now available high-quality data from LEAR and KEK, 
we have been able to extend the methods used in these AW PWAs to perform an energy- 
dependent PWA of all pp scattering data below p^ = 925 MeV/c (T^ = 379 MeV). This 
work was started in 1987 P] and has only recently been finished 0. The same methods of 
PWA have also been applied ||^ to the strangeness-exchange reaction pp AA, for which 
the PS185 group at LEAR has obtained beautiful data. In the next section we review the 
theoretical ideas behind these Nijmegen PWAs, and in section III we apply these ideas 
and methods to the case of pp scattering. In section IV some results of this pp PWA are 
presented and discussed. 



After almost a decade of LEAR, it is fair to say that in this field theory has some 
catching up to do with respect to experiment. Since the partial- wave amplitudes or the 
phase-shift parameters are in a sense the meeting ground between theory and experiment, 
the results of the present PWA should be very useful in many ways. They can be used 
to improve models |l^ for the NN interaction. Apart from the fact that this provides 
independent and complementary 0] information about the spin- and isospin structure of 
the NN force, the NN interaction is needed as input in many other pp subfields. Studies 
of for instance protonium (the pp atom) or specific annihilation processes like pp —>■ tt+tt^, 
K^K~ require a realistic treatment of the initial pp interaction. At the same time, this 
PWA could be helpful in planning new experiments at LEAR, the future of which is of 
course crucial to this field. 

II. Methods of partial-wave analysis 

The hallmark of the Nijmegen energy-dependent PWAs is the sophisticated manner in 
which the energy dependence of the partial- wave amplitudes is parametrized. At the basis 
of the PWA is the trivial observation that in the low-energy region (long wave lengths) 
the long-range interaction is very important. It is this long-range interaction that is 
responsible for the rapid variations with energy of the scattering amplitudes. Short-range 
interactions lead to much slower energy variations of the amplitudes. One usually looks 
for a function in the problem that one can parametrize as easily as possible, i.e. one 
that does not contain the contributions from these long-range processes. Because these 
long-range interactions are at the same time model independent (in the sense that they 
are or at least should be the same in all NN and AW models), they can then be taken 
into account separately and exactly. 

It is, of course, not a good idea to try to parametrize the partial-wave 5* matrix 
itself, since it does contain all of these long-range effects. As a function of complex 
energy, the S matrix has a (kinematical) right-hand unitarity cut, other right-hand cuts 
due to the coupling to inelastic channels, and (dynamical) left-hand cuts due to particle 
exchanges. The left-hand cuts that are the closest to the origin = correspond to 
the longest-range processes. The left-hand cuts that start far away from the origin are 
due to the short-range interactions. For instance, the infinite-range Coulomb potential 
iy ~ l/r) produces an essential singularity and a branch point at = 0, vacuum 
polarization {V ~ exp(— 2mer)/r^/^) produces a cut at = —0.6 keV, and one-pion 
exchange (V ~ exp{—mT,r)/r) leads to a cut starting at Tl = —9.7 MeV. One sees that 
the crux is to find a quantity in which the cuts nearest to the origin are not present. This 
quantity then allows an analytical parametrization in energy or k"^ in an enlarged domain 
up to the next left- or right-hand cut present. 

A familiar example of such a quantity with improved analyticity properties is the 
modified effective-range function |ll|. The Coulomb-modified effective-range function 
for the pp ^Sq state was originally derived (in a rather intuitive way) by Landau and 
Smorodinsky [|12]. When only the Coulomb potential is present the boundary condition 



for the radial wave function $(0) =0 is of course satisfied by F, the regular Coulomb 
wave function (for i = 0). Suppose that there is an additional short-range interaction. 
When the wave length is very large (very low energy), one can take the limit in which 
the range of this additional (strong) interaction goes to zero. Then its presence is only 
revealed by a modified boundary condition at r = 0, which is now satisfied by a linear 
combination of F and G, the irregular Coulomb wave function 

$(r) = F{r) cot 6o + G{r) , (1) 



where So is the nuclear phase shift in the presence of the Coulomb interaction ((5(^5*0) = 
So + ctq), as can be seen from the asymptotic behavior of $. An equation for cot^o can 



then be obtained by evaluating the logarithmic derivative of the wave function, which 
we call P{k,r), for A; — > 0. In the np case this quantity P{k,0) = fccot^o approaches a 
constant: 



In the pp case, the evaluation has to be done at r = e because of a term In e that appears 
due to the singular behavior of G. This term one absorbes in the constant — 1/a, along 
with some further constant terms. Then one lets e — > and immediately obtains the 
Coulomb- modified ("zero-range") effective-range function. It can be shown that after 
these manipulations the resulting left-hand side of Eqn. (^) is an analytical (actually 
meromorphic) function of the energy, so that the right-hand side can be written as a 
power series in k"^ (this means dropping the zero-range approximation). 

The analytical expansion of the Coulomb-modified effective-range function breaks 
down already at = —9.7 MeV, where the one-pion-exchange cut starts. It is pos- 
sible to derive a new "pion-modified" effective-range function from which also this cut 
has been removed 0]. Let the regular and irregular wave functions for the case where 
only the Coulomb and pion-exchange potentials are present be called and G^. (For 
the purpose of the present discussion, we ignore vacuum polarization.) The wave function 
can then be written as 



where 6^ is now the phase shift due to the short-range remainder of the strong interaction 
((5(^5*0) = 6o + ttq + (To, where ttq is the one-pion-exchange phase shift in the presence of 
the Coulomb potential). However, proceeding in similar fashion as above, one encounters 
an important problem here. The evaluation of P{k,e) has to be done numerically, since 

and are not known in analytical form. Due to the singular behavior of when 
£ ^ 0, it is very hard to maintain sufficient numerical accurary, especially for higher 
orbital angular momenta. 

At this point one has to realize that this numerical problem of the modified effective- 
range function is really an artificial problem: it crops up due to the singular behavior of 
the irregular function of the long-range potential near the origin. However, it is precisely 
this short-range interaction that one wants to parametrize, since it is essentially unknown, 
very complicated, and leads to only slow energy variations of the scattering amplitudes. 
Looking at Eqn. (^, one observes that it is valid for any r, so why not evaluate P{k,r) 
at a finite value r = b, instead of at r = e? 

This is essentially what is done in the Nijmegen PWAs. The wave functions are ob- 
tained by solving the (relativistic) Schrodinger equation. Suppose one starts at a point Too 
where only the Coulomb potential is present. Integrating inwards, one picks up sequen- 
tially the contributions (varying rapidly with energy) from the electromagnetic potentials, 
one-pion exchange, and contributions (varying slower with energy) from other meson ex- 
changes. One then stops at a point r = b. If there are no additional interactions for 
r < b, the boundary condition P{k, 6) at r = 6 is obviously satisfied by the regular wave 
function corresponding to precisely this potential tail. For small enough b the model used 
for r > b will of course not be correct, and the boundary condition has to be modified, 
as in the above examples. In practice, it works the other way: one starts integrating 
at r = 6 and P{k,b) is parametrized as a function of energy. Also the best value for b 
is determined by fitting the data. In general multichannel problems P{k, b) becomes a 
matrix. This P matrix has the required improved analyticity properties. When there 
are no nearby right-hand cuts, it is an analytical (again: actually meromorphic) function 
of fc^ in a domain bounded by the nearest left-hand cut not removed by including (or 
including incorrectly) the corresponding exchange in the potential tail for r > b. It can 




(2) 




(3) 



happen, of course, that short-range dynamics gives rise to a rapid energy variation of the 
amphtudes, as in the case of a resonance. This would have to be taken into account in 
the P matrix, for instance by including a pole in the parametrization. It is seen that the 
formalism used in the Nijmegen PWAs is similar to the boundary-condition approach to 
the strong interactions that goes back to the work of Feshbach and Lomon [|l^ and earlier. 
The philosophy, however, is very different. The term P matrix (for "pole" matrix) was 



introduced by Jaffe and Low |]T4[ in the framework of the bag model. 



III. An antiproton-proton partial-wave analysis 

Let us now be more specific and apply the foregoing ideas to the case of pp scattering. In 
all the Nijmegen PWAs, the two-body scattering process is described with the relativis- 



tic Schrodinger equation ||I5|, which is essentially a coordinate-space version of the 



Blankenbecler-Sugar equation. It reads the same as the ordinary Schrodinger equation 

(a + _ 2mV) ^{r) = , (4) 

except that the proper relativistic relation between energy and momentum is used. It is 
well known how to derive the potentials for use in this equation ||T^, IB]. In this relativistic 



framework, there is no known quantum-mechanical interpretation for the "wave function" 
ip{v). It is perhaps best to regard it as just a tool that allows one to compute the correct 
relativistic scattering amplitude (e.g. the poles are the correct bound states). We solve 
Eqn. (^) for the coupled pp and Tin channels. The mass difference between proton and 
neutron is included in order to account for the fin threshold at pl = 99 MeV/c. 

The interaction in the region r > 6 is described by a theoretically well-founded NN 
potential. This potential is given by 

V = Vc + Vmm + Vn , (5) 

where Vc and Vmm are the relativistic Coulomb and magnetic-moment interaction re- 
spectively. Vn is the NN meson-exchange potential. It consists of one-pion exchange and 
the (charge-conjugated) heavy-meson and pomeron exchanges from the 1978 Nijmegen 
NN potential [0. As argued in the previous section, the rapid energy variations of the 
amplitudes due to the long-range electromagnetic interactions and one-pion exchange are 
now included exactly. 

Let us next turn to the parametrization of the short-range interactions for r < 6 by 
way of the P-matrix boundary condition at r = 6 = L3 fm. Due to the coupling to the 
annihilation channels, the S matrix has a right-hand cut starting already to the left of 

= 0. (In the pp case this cut starts only at the pp ppir^ threshold at = 280 
MeV.) As these annihilation processes are of short range (and so give rise to slow energy 
variations of the amplitudes), this right-hand cut has to be present in the P matrix, which 
we therefore take to be complex. (Similarly, the effective-range parameters for the NN 
case are complex.) The choice of the value for b is rather critical, more so than in the NN 
case (where it was taken to be 6 = 1.4 fm). The best results are obtained for 6 = 1.3 fm. 
Since for r > 6 we use only a real potential, the coupling to the annihilation channels is 
completely represented by the boundary condition. We conclude therefore that the range 
of the annihilation process is in fact about 1.3 fm [0. 

The electromagnetic interactions that we use are adapted from the improved Coulomb 
potential ||16[. This potential, designed specifically for use in the relativistic Schrodinger 



equation, contains relativistic corrections to the static Coulomb potential and (in its off- 
shell behavior) the main contributions from the two-photon-exchange diagrams. All these 
effects are included in the Nijmegen pp PWA @, |^, as well as the vacuum-polarization 



potential. In our case it suffices to use the following spin-dependent one-photon-exchange 
potentials 
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The magnetic moments of the proton and neutron are fip = 1 + Kp = 2.793 and /i„ = 
Kn = —1.913, respectively. The use of a' in the central potential for pp pp takes care of 
the main relativistic corrections to the Coulomb potential. It is given by a' /a = 2k/MvL 
where vl is the velocity of the antiproton in the laboratory system. At 600 MeV/c one 
has for instance vl = 0.54 and a' /a = 1.135, a correction of 13.5% to the static Coulomb 
potential. The spin-orbit potential comes from the interaction of the magnetic moment 
of one particle with the Coulomb field of the other particle (and is consequently absent in 
nn —>■ nn) . It includes a relativistic correction due to the Thomas precession. The tensor 
potential comes from the interaction of the two magnetic moments. Vacuum polarization 
and two-photon-exchange effects are negligible in our case. The proper treatment of 
these electromagnetic effects in the evaluation of the scattering amplitudes is a nontrivial 
matter [0. The following simple one-pion-exchange potential without a form factor is 
used 
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where m is the mass of the pion and f^Nn — 0.0745 is the rationalized pion-nucleon 



coupling constant 18 



The mass difference between the 7r° and vr^ is included. 
Let us finish this section with some more general remarks about PWAs. Even for 
the pp case, where the database is of high quality and the observables are very well 
mapped out, a PWA is impossible without a substantial amount of theoretical input or 
constraints. For the np and pp PWAs, this is true a fortiori. For instance, one has to 
make some assumptions about the validity of symmetries like charge independence or (as 
in our case) charge conjugation. Obviously, one has to careful here: sometimes general 
physical principles are inspired by local renormalizable field theories and not strictly valid 
for extended objects like hadrons. A good example can be found in ttN PWAs, where one 
usually implements full Mandelstam analyticity |jT9[. The amplitudes are assumed to be 
analytic functions of the two complex variables s and t except for singularities from the 
mass spectrum and unitarity. These amplitudes then exhibit crossing symmetry. It is not 
clear at all to what extent low-energy hadron dynamics actually satisfies this symmetry. 

Using strong and mostly model-independent theoretical constraints it has turned out to 
be possible to perform an energy- dependent or multienergy PWA of the pp data. However, 
it is quite a different ballgame to perform energy-independent 01 single- energy pp PWAs. 
In a single-energy pp PWA one has to determine in principle 20 phase-shift parameters for 
each J 7^ (8 for J = 0), which is four times as many as in a single-energy np PWA 0! 
Almost certainly the present pp database does not allow satisfactory energy-independent 
PWAs. One has to realize, however, that even in the NN field the usefulness of energy- 
independent PWAs is more limited than is perhaps generally thought. When one has a 
(?ooc/ energy-dependent PWA, the best value for a phase shift (or the pion-nucleon coupling 
constant!) is definitely the one determined in the energy- dependent PWA, and not the 
one from an energy-independent PWA. One reason is that an energy-independent PWA 
contains no information about the energy dependence of the amplitudes. This makes it for 



instance less stable than an energy-dependent PWA with respect to the addition of new 
data to the database. Also, a set of energy-independent PWAs is usually overparametrized 
compared to a good energy- dependent PWA in the same energy region. It thus almost 
certainly contains noise. For an extensive discussion of this important point, see Ref. 

IV. Some results of the analysis 

While in AW PWAs there is essentially agreement on the correct database (especially for 
pp), we had to spend a lot of time and effort into collecting, scrutinizing, and cleaning 
up the world set of pp scattering data, which contains a lot of flaws and contradictory 
data. Exactly the same statistical arguments were used in this process as were used in the 
set-up of the Nijmegen NN database 0, ^. This means for instance that data with a very 
improbable high or low are rejected on statistical grounds. The resulting Nijmegen 
1993 pp database in the momentum interval 119-923 MeV/c is unique in the world and 
consists of Wdata = 3646 pp data. It is extensively discussed in Ref. |^. In the final fit 
to this database we reach = 3801.0 or x^/^data = 1.04. The number of boundary- 
condition parameters needed is 30, which is a reasonable number, in view of the fact that 
21 parameters were needed in the Nijmegen pp PWA and an additional 18 in the np PWA. 
The total number of degrees of freedom is A^^df = 3503, which means that x^/^df = 1.09. 

If the database is a correct statistical ensemble and if the theoretical model is correct, 
one expects that (x^) = A'df = 3503 with an error of ^/2Nd{ = 84. This means that in 
our PWA we are 298 or only 3.5 standard deviations away from the expectation value for 
X^. We conclude that although there is still room for improvement, our energy- dependent 
solution is essentially correct statistically. As a consequence, the values for the phase- 
shift parameters (and also for the pion-nucleon coupling constant) and the statistical 
errors (obtained in the standard manner from the error matrix) are essentially correct. 

In our 1991 preliminary PWA we were able to determine the charged-pion-nucleon 
coupling constant = fpnn+ fnpw- from the data on the charge-exchange reaction 
pp nn, in which only isovector mesons can be exchanged. The result found was 

= 0.0751(17), at the pion pole. The error is purely statistical. In our final analysis, we 
have repeated the determination of but this time from the complete 1993 Nijmegen 
database. The coupling constants of the neutral pion were kept at the value of fpp^o = 
/2^^o = 0.0745 ig. We now find = 0.0732(11), at the pion pole. This resuh supersedes 
our previous value from Ref. . Again, the error is of statistical origin only. In view of the 
enormous amount of work involved, it is very hard to estimate possible systematic errors 
on this result. We have checked that there are no systematic errors due to form-factor 
effects or due to uncertainties in p^(770) exchange. In the Nijmegen pp PWA systematic 



errors could be more thoroughly investigated and they were found to be small |jT8|. In 
our case the systematic errors are probably larger than for the pp case, but it is very 
encouraging that the result for is in good agreement with recent determinations = 
0.0735(15) from n^p ^ scattering and o = 0.0745(6) and = 0.0748(3) from NN 



scattering |TB[. Very probably the new LEAR experiment PS206 on pp nn will further 
constrain the charged-pion-nucleon coupling constant. 

In Fig. |I| the differential cross section at 693 MeV/c and the analyzing power at 
875 MeV/c are shown for pp nn. The data are from PS199 0. One can see the truly 
remarkable accuracy of the cross-section data and the analyzing-power data in the forward 
region. The "dip-bump" structure in da/dQ at forward angles is due to the interference 
of one-pion exchange and a smooth background. 

The fact that the available charge-exchange data already pin down the charged-pion 
coupling constant with a remarkable small statistical error is only one example of how at 
present quantitative information can be extracted from the pp system. We can mention 




Figure 1: Differential cross section at 693 MeV/c and analyzing power at 875 MeV/c for 
the charge-exchange reaction pp — > Tin. The data are from PS199 [0. The curves are 
from the Nijmegen PWA [0. 



some subtle effects that are also visible in the data. The use of a' instead of a, i.e. the 
main relativistic correction to the static Coulomb potential, gives a drop of A^^ = 30, or 
5.5 standard deviations. The inclusion of the magnetic-moment interaction gives a drop 
of Ax^ = 14, or 3.7 standard deviations. Even the use of the correct pion masses instead 
of an average mass of 138 MeV is a 3 standard-deviation effect. 

Since the present pp PWA is the first of its kind, we have also proposed a convention 
for extracting phase-shift and inelasticity parameters from the S matrix. In the presence 
of coupling to annihilation channels the S matrix describing AW scattering is only a 
submatrix of the much larger multichannel 5* matrix. It is therefore still symmetric, but 
no longer unitary. This doubles the number of parameters needed. For the partial waves 
with £ = J, s = 0, 1 one obviously writes S = riexp{2i6). For the states with i = J ± 1, 
s = 1, coupled by the tensor force, six parameters are needed to parametrize the 2x2 
S matrix. In this case it is not so easy to think of a convenient parametrization which 
satisfies all constraints from unitarity, is completely general, and free from nontrivial 
ambiguities. We have used a generalization of the "bar-phase" convention commonly 
used in NN scattering. One writes (with the notation 6ej for the phase shift) 



cos Ej ismej 
i sine J cos e j 



cos Ej I sm sj 
i sin Ej cos e j 



exp(i(5) 



(9) 



where 6 = diag((5j_i,j, ^ j+i,j) and ej is the mixing parameter. H'^ is a three-parameter 
real and symmetric matrix representing inelasticity. It can be diagonalized in Blatt- 
Biedenharn fashion 



cos ujj — sm ujj 
sin uj cos Uj 



VJ-i,J 
r]J+i,J 



cos ujj sm u J 
- sin cuj cos u J 



(10) 



where the diagonal matrix contains the "eigeninelasticities" rjj^i^j and rjj+i^j, and uj is 
again a mixing parameter. We are presently in the process of doing a careful evaluation of 
these phase-shift and inelasticity parameters and their errors. Unfortunately, this involves 
a large amount of work. These and other issues will be the subject of future publications. 
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